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FOREWORD 


1.  This  is  a report  documenting  the  demilitarizing  of  an  HBX-1 
loaded  Mk  25  Mod  1 Mine  by  an  in  situ  autoclave  melting  procedure 
and  an  analysis  of  why  the  technique  was  aborted  and  the  thermo- 
chemistry involved. 

2.  The  effort  reported  herein  was  authorized  and  funded  under  the 
Naval  Sea  Systems  Command  (SEA-9912)  Work  Request  63601  of  6 Nov  1975. 
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THERMOCHEMISTRY  AND  THE  DEMILITARIZATION  OF  EXPLOSIVES 


I.  INTRODUCTION 

Prior  to  the  national  awareness  of  a limited  energy  supply,  the 
main  thrust  of  governmental  operations  concerned  with  the  manufacture, 
use  and  disposal  of  explosive  ammunition  was  to  meet  or  better  federal, 
state  and  local  pollution  and  environmental  standards.1  2 Rigid 
restraints  were  placed  on  deep-sea  dumping  and  open  air  burning  as 
acceptable  means  for  disposal  of  obsolete  or  substandard  ammunition. 

As  a result,  interest  was  revived  in  procedures  for  reclamation  and 
recovery  of  explosive  moieties  for  either  reuse  or  commercial  sale. 

The  majority  of  these  processes  were  based  upon  either  steam  or  water 
washout  of  the  explosive  components  which  were  then  recovered  by  a 
variety  of  decantation  and  filtration  techniques.  Without  recycling 
or  reuse  of  the  water  effluents,  these  processes  were  both  costly  and 
energy  consumptive. 

The  oil  embargo  of  the  early  1 970 ’ s brought  the  need  for  energy  con- 
servation into  sharp  focus.  The  U.S.  Navy's  demilitarization  programs 
were  expanded  to  emphasize  energy  conservation  as  well  as  pollution 
control.  In  keeping  with  this  energy  conservation  mandate,  the  Naval 
Explosives  Development  Engineering  Department  (NEDED),  of  the  Naval 
Weapons  Station  (NAVWPNSTA),  Yorktown  demonstrated  on  a pilot  scale 
that  meltable,  TNT-based  explosives  could  be  recovered  safely  by  wash- 
out with  recycled,  pressurized  hot  water.* * 3 

As  a further  extension  of  the  energy  conservation  concept.  Battel le 
Columbus  Laboratories,  working  under  contract4  to  the  Western  Division 
of  the  Naval  Facilities  Engineering  Command,  had  calculated  that  even 
further  energy  conservation  could  be  realized  by  completely  melting 
TNT-based  explosive  in  situ  in  the  respective  ammunition.  Battelle's 
concept  was  to  place  the  ordnance  item  in  an  autoclave  heated  by 


National  Environmental  Policy  Act  (PL-91 -90). 

Executive  Orders  11507  & 11514. 

3NAVWPNSTA  Yorktown  files  (NEDED  File  #897  Uun  1971 -Sep  1972). 

4Kim,  B.C.,  et  al .,  Battelle  Columbus  Laboratories,  Final  Report  on 
Explosive  Melt  Out  Study  (N62474-74-C3913),  Sep  1975. 
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10  to  15  pounds  per  square  inch  gage  (psig ) recycled  steam  until 
melting  was  completed.  The  molten  explosive  was  then  to  be 
syphoned  from  the  warhead,  and  either  delivered  directly  to  a 
solidifying  flaker  belt  and,  hopefully,  be  packageable  for  reuse 
and/or  sale  without  further  processing  or  else  delivered  to  an 
incorporation  kettle  for  reconstitution  to  specification  grade 
standards. 

For  example,  from  preliminary  experiments  at  the  Naval  Ammunition 
Depot  (NAD),  Hawthorne  the  energy  requirements  for  autoclave  demili- 
tarizing TNT  and  HBX-1  loaded  Mk  25  Mod  1 Mines  were  estimated  to  be 
an  order  of  magnitude  (10X)  lower  than  the  energy  requirements  for 
either  a steam-out  or  washout  process.4* 

Prior  to  the  NAD  Hawthorne  experiments,  Allegheny  Ballistic  Laboratory 
(ABL)  had  been  contracted  to  run  a hazard  and  risk  analysis  on  the 
proposed  process.  Their  calculations  and  estimates  were  that  the 
probability  of  an  explosion  using  the  proposed  autoclaving  procedure 
for  specified  time  periods  were  less  than  5*10-9.5  It  was  on  the 
basis  of  these  estimates  that  the  Hawthorne  experiments  were  conducted. 

NEDED  is  chartered  to  either  develop  or  approve  all  processes  relevant 
to  the  loading  or  demilitarization  of  main  charge  explosives  in  Navy 
ammunition  by  Naval  Sea  Systems  Command  instruction  and  was  directed  to 
evaluate  the  autoclaving  process  in  more  detail.  This  report  concerns 
itself  with  that  evaluation,  and,  in  addition,  summarizes  some  signifi- 
cant data  from  the  thermochemistry  of  explosives  literature.  The  prime 
objective  is  to  impress  explosive  process  and  design  engineers,  not 
normally  conversant  with  the  thermochemistry  of  explosives,  with  the 
limitations  as  well  as  the  value  of  applied  thermochemical  explosion 
theory. 

II.  EXPERIMENTAL  OBJECTIVES  OF  THE  AUTOCLAVE  MELTING  PROCESS 

NEDED 's  objectives  essentially  were  to  verify  and  refine  the 
findings  at  NAD  Hawthorne  listed  in  the  Battel le4  report.  Namely  that: 

• Gross  steam  consumption  for  autoclave  melting  explosives 
in  Mk  25  Mod  1 TNT  and  HBX-1  loaded  Mines  ranged  from  0.24 

to  0.30  pounds  of  steam  per  pound  of  explosive,  respectively. 

• Net  steam  consumption  could  be  reduced  to  as  low  as  0.12 
pounds  of  steam  per  pound  of  explosive  by  reducing  heat 
losses  from  the  melt-out  chamber  by  better  insulation. 


kIbid.  (*p  14). 

5Sawyer,  S.O.,  & Williams,  G.M.  Allegheny  Ballistics  Laboratory 
Final  Rept  No.  A0827-520-03-006,  Hazards  Analysis  Report,  Part  II, 
Jan  1976  (prep  for  Battel le  Columbus  Laboratories). 
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• Melting  times,  using  10  to  15  psig  steam  and  without 
agitation,  were  somewhere  between  6 and  13  hours  for  the 
two  mines  respectively  or,  in  any  case,  less  than  16  hours. 
This  latter  objective  was  considered  important  because  of 
the  previously  referenced  ABL  hazards  estimates.5 

• Molten  hot  melt  could  effectively  be  skimmed  off  from  the 
bulk  of  the  remaining  explosive. 

• Molten  explosive  could  be  syphoned  directly  from  the  war- 
head for  immediate  pelletization  or  for  reconstitution  if 
necessary. 

The  experimental  setup  at  the  NAVWPNSTA  Yorktown  was  essentially 
equivalent  to  that  at  NAD  Hawthorne  as  described  in  Battel le's  report 
and  as  illustrated  in  Figures  1,  2 and  3. 

A Mk  25  Mod  1 HBX-1  loaded  Mine,  containing  approximately  1200  pounds 
of  explosive  was  sealed  into  the  autoclave  and  the  following  sequence 
of  events  noted  and  carried  out: 

• At  0000  hours,  15  psig  steam  was  applied  to  the  autoclave 
jacket. 

• At  0200  hours,  approximately  20  pounds  of  hot  melt  had  exuded 
up  into  the  vent  tube  and  was  removed  with  a dip  stick. 

• At  0830  hours,  probing  via  the  vent  tube  showed  that  explosive 
was  melted  through  to  the  bottom  of  the  mine  but  that  a mass 
of  explosive  against  the  "lug"  side  of  the  mine  still  remained 
solid.  Some  bubbling  or  gassing  was  observed  at  this  time. 

• At  1130  hours,  melting  was  not  complete  and,  despite  inter- 
mittent stirring,  gassing  had  increased  to  a steady  and 
excessive  degree.  Accordingly,  the  decision  was  made  to 
abort  the  experiment.  The  steam  to  the  autoclave  was  cut  off 
at  this  time  and  the  latter  vented. 

• At  1215  hours,  approximately  250  pounds  of  molten  explosive 
was  syphoned  into  a reprocessing  kettle.  Steady  bubbling 
continued  in  the  molten  explosive  remaining  in  the  mine. 
Therefore,  cooling  of  this  residual  material  was  hastened 
by  removing  the  "top-hat"  and  lid  to  the  autoclave  and  by 
dumping  approximately  600  additional  pounds  of  the  molten 
explosive  into  pellet  trays. 


5Ibid. 
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FIGURE  3.  AUTOCLAVE  CLOSED  READY  FOR  HEATING 
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• No  analysis  was  made  of  the  escaping  gas  but  it  was  the 
certain  conviction  of  the  cognizant  engineer  that  the  gas 
was  a reaction  product  and  not  the  release  of  entrapped  air. 

• Laboratory  analyses  of  the  syphoned  HBX-1  (Sample  No.  1), 
HBX-1  dumped  into  pellet  trays  (Sample  No.  2),  and  a 
standard  HBX-1  formulation  (Sample  No.  3)  are  shown  in 
Table  I. 


TABLE  I.  MK  25  MOD  1 MINE  HBX-1  ANALYSES 


Sample 

No. 

Composition  (%) 

Impact 

sens iti vitya 
(cm) 

VTSb 
(ml /g ) 

H20 

TNT 

Hot 
mel  t 

A1 

Wax 

RDX 

CaCl2 

1 

0.2 

40.8 

0.3 

14.9 

2.6 

40.8 

0.4 

84.4 

1.0 

2 

0.2 

37.7 

0.4 

17.1 

4.0 

40.7 

0.5 

82.2 

0.9 

3 

0.2 

37.8 

0.0 

17.1 

4.7 

39.7 

0.5 

85.9 

0.1 

a50%  ht:  Bruceton  Method/NOL  Machine/Type  12  tools/2.5  kg  wt/ 

25  trials/35  +_  2 mg/  5/0  sandpaper. 

^Vacuum  Thermal  Stability:  ml/g/48  hrs/100°C/STP. 


A highly  significant  point  to  be  noted  from  the  vacuum  stability 
results  is  the  tenfold  increase  in  gas  evolved  from  remelted  material 
versus  the  virgin  standard.  Further,  methyl  violet  tests  of  these 
vacuum  stability  gas  products  from  the  reclaimed  HBX-1  were  positive 
for  nitrogen  oxides,  whereas  tests  were  negative  for  newly  batched 
material  tested  under  the  same  conditions.  These  factors  will  be 
readdressed  later. 

Consequent  to  this  experiment,  a thorough  review  of  the  thermochemistry 
of  explosives  literature  was  made.  Section  III  summarizes  the  develop- 
ment of  the  thermochemical  theory  of  explosions.  Section  IV  reviews 
its  applications  and  limitations. 
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III.  THERMOCHEMISTRY  OF  EXPLOSIVES 

Some  50  years  ago,  N.  N.  Semenov  promulgated  the  theory  of 
thermal  explosion.  Since  that  time,  steady  progress  has  been  made  in 
developing  mathematical  tools  with  which  both  explosive  process  and 
weapon  design  engineers  are  able  to  calculate  "critical  temperatures" 
(Tc)  and  "times-to-explosion"  (te)  for  explosive  compositions  having 
specific  geometric  configurations  and  temperature  environments. 

Differences  between  calculated  and  experimentally  determined  values 
for  Tc  and  te  for  "chemically  pure"  explosive  systems  have  been 
steadily  narrowed.  This  progress  is  due  primarily  to  two  factors: 

1)  more  exact  solutions  to  the  heat  conduction  equation  for  thermal 
explosion  have  been  developed,  and  2)  more  accurate  kinetic  constants 
which  are  used  in  the  heat  balance  equation  have  been  generated  for 
specific  explosive  systems. 

Several  excellent  reviews  relevant  to  this  progress  have  been  pub- 
lished and  it  is  not  the  purpose  here  to  add  to  that  list.  Rather 
our  objective  is  to  show  that  as  agreement  between  theory  and  ex- 
periment is  approached  for  chemically  pure  explosive  systems,  wider 
recognition  should  be  given  to  major  discrepancies  between  calculated 
and  observed  Tc  and  te  values  (particularly  the  latter)  obtained  from 
production  grade,  adulterated,  or  less  pure,  explosive  systems. 

While  this  note  of  caution  will  come  as  no  surprise  to  long  time 
investigators  in  this  field,  it  is  being  reemphasized  here  for  two 
reasons.  The  first  is  so  that  explosive  engineers  realize  that  in 
using  the  kinetic  constants  for  pure  explosives  and  in  making  simpli- 
fied assumptions  about  the  type  of  reaction  mechanism  occurring,  they 
are  calculating  the  "most  optimistic"  times-to-explosion.  The  second 
reason  is  to  encourage  an  expanded  and  accelerated  effort  in  the 
research  and  development  community  to  further  identify  and  treat  com- 
peting reaction  mechanisms  that  will  lead  to  formulae  from  which 
"worst  case"  predictions  can  be  made  with  higher  confidence. 

Therefore,  the  theory  of  chemical  explosion  is  resurrected  here  only 
to  the  degree  felt  necessary  for  the  uninitiated  explosive  engineer 
to  gain  some  insight  into  its  derivation  and  application.  The 
summary  of  methods  for  calculating  critical  temperatures  and  times- 
to-explosion  has  been  abstracted  from  several  sources.  Key  among 
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them  are  those  now  referenced.6  7 8 9 10  11  Pioneering  efforts  by 
Frank-Kamenetskii , Chambre,  Carslaw,  Jaeger  and  Eyring,  to  name  but 
a few,  already  have  been  liberally  cross-referenced  among  the  publi- 
cations cited  and  are  not  listed  as  prime  references.  For  those 
engineers  less  interested  in  thermal  explosion  theory  than  in  how 
and  why  it  may  vary  from  experimental  data,  it  is  recommended  they 
skip  to  Section  IV. 

A.  Thermal  Explosion 

Materials  that  undergo  self-heating  by  exothermic  decompo- 
sition will  do  so  at  accelerated  rates  as  the  temperature  rises.  If 
conduction  at  the  material  boundary  is  not  sufficient  to  dissipate 
the  heat,  thermal  explosion  will  eventually  occur.  The  highest 
boundary  or  surface  temperature  at  which  a "steady  state"  can  exist 
is  that  wherein  heat  losses  and  heat  production  are  in  equilibrium. 

This  temperature  is  defined  as  the  critical  temperature,  i.e.,  the 
temperature  which  if  exceeded,  will  result  in  a thermally  induced 
explosion.  Obviously,  since  heat  flow  is  boundary  controlled,  critical 
temperature  is  strongly  influenced  by  the  geometry  and  size  of  the 
reacting  entity. 

This  phenomenon  of  self-heating  by  many  materials,  including  explosives, 
can  be  described  by  the  heat  conduction  equation  (neglecting  convective 
heat  flow)  as: 

Self-heating  = Heat  Loss  by  Conduction  + 

Production  of  Heat  by  Chemical  Reaction 

Quantitatively,  this  expression  is 


pc 


3T 

9t 


av2t  + PQZirE/RT 


0) 


6Longwell,  P.A. , NAVWEPS  Rept  7646,  Calculation  of  Critical  Tempera- 
ture and  Time  to  Explosion  for  Propellants  and  Explosives,  Mar  1 961 . 
7AMCP  706-180,  Principles  of  Explosive  Behavior,  Apr  1972. 

8Zinn , J.  & Mader,  C.L.,  J.  Appl . Phys.  H (1960)  323. 

9Rogers,  R.N.,  Thermochimica  Acta  XL  (1975)  131-9. 

10Hardee,  H.C.,  Lee,  D.O  & Donaldson,  A.B.,  Combustion  and  Flame 
]8  (1972)  403-10. 

nEnig,  J.W. , Proc.  Roy.  Soc.  305A  (1968)  205. 
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where  p = density  g/cc 

c = specific  heat  cal/g°K 
T = temperature  prior  to  explosion  °K 
A = thermal  conductivity  cal/cm°Ksec 

v * Laplacian  operator  in  Cartesian  space  coordinates  x,  y, 

i.e.  = »!_  + ii  + iL 
9x2  ay2  az2 


z. 


t = time  sec 

Q = heat  of  reaction  per  unit  mass  cal/g 
Z = pre-exponential  factor  sec-1 
E = activation  energy  cal /mol 
R = gas  constant  1.987  cal/mol°K 


Proximate  and  exact  solutions  to  Eq.  (1)  have  been  derived  with 
increasing  ingenuity  since  its  first  proximate  solution  by  Frank- 
Kamenetskii  for  the  steady  state  condition  (i.e.,  where  ^-  = 0)  of  a 
slab  with  both  faces  held  at  constant,  adiabatic  temperature. 


In  general,  for  symmetric  geometries,  conduction  depends 

single  space  coordinate  so  that  v2  reduces  to  — + — — 1 

x2  x\ax| 

1 or  2 for  planar,  cylindrical  or  spherical  symmetry  and 
controlling  space  coordinate;  e.g.,  for  a sphere  x = the 

j\T 

and  m = 2 and,  where  9c  =0,  Eq.  (1)  reduces  to 


only  on  a 
where  m = 0,  • 
x is  the 
radius,  r, 


pQZfE/RT 


0 


or,  more  commonly,  letting  0 = 


d2e  + m/del  + pQZR^-i/q  _ Q 
dx2  x d x]  AE 


(2) 


(3) 
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Until  recently,  exact  solutions  to  the  nonlinear  Eq.  (3)  had  not 
been  developed.  However,  proximate  solutions  to  Eq.  (3)  have  been 
derived  by: 


• making  the  assumption  that  the  temperature  of  the  explosive, 
T,  just  prior  to  explosion  is  very  close  to  a fixed  reference 
temperature,  T0,  such  as  the  boundary  temperature,  or  in  a 
specific  case,  the  critical  temperature. 


i .e. 


L,±i 

' RT  RTol 


where  L_I°  <<  1 
To  / T0 


• defining  a new  temperature  variable 

0 = -^  (T-T o ) and  de  = -^  dt 

Rll  rjI 

• defining  a dimensional  positioned  coordinate,  c,  as  ? = x/a 
where  a is  the  significant  dimension  of  the  explosive,  such 
as  half  of  the  thickness  of  a slab  or  the  radius  of  a 
cylinder  or  a sphere. 


Appropriate  substitutions  of  these  values  in  Eq.  (3)  result  in 


d2  0 
dt2 


(4) 


This  proximate  equation  has  been  solved  for  various  geometries  and 
boundary  conditions.  Of  particular  importance  is  the  parameter  in 
the  above  equation  generally  identified  as  6,  the  so-called  "shape 
factor,"  where 


5 = i!^£2|it-E/RT0 
xrtJ 


(5) 


(At  the  critical  temperature,  T0  = Tc.) 
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The  maximum  or  "critical"  value  for  the  parameter  6,  i.e.,  6cr,  was 
calculated  by  early  investigators  to  be  a constant  dependent  upon 
size  and  shape;  i.e.,  Scr  = 0.88,  2.00  and  3.32  for  the  infinite  slab, 
cylinder  and  sphere,  respectively.  Zinn  and  Mader  solved  Eq.  (4) 
numerically  for  these  geometries  with  nonadiabatic,  nonsteady  state 
systems.  Later,  Enig,  Longwell,  Conrad,  Hardee,  et  al.%  calculated 
6 values  for  other  boundary  conditions  for  the  significant  dimension 
of  the  explosive  body  and  showed  that  the  6cr  values  just  listed  are 
limiting  values.  Ultimately,  Shouman,  Donaldson  and  Tsao12  offered 
an  exact  solution  for  the  heat  transfer  equation.  However,  since 
many  of  the  sizes  and  shapes  of  high  explosives  systems  of  interest 
to  the  design  engineer  have  6 values  that  approach  those  values 
listed,  their  utility  is  of  general  application  and  therefore,  Eq. 

(5),  as  we  shall  next  see,  is  valid  for  use  with  simple  systems. 

B.  Critical  Temperatures 

For  the  cases  just  described,  where  6cr  - 0.88,  2.00  or 
3.32,  Eq.  (5)  [which  in  fact  is  the  integrated  solution  to  Eq.  (4) 
for  slabs,  cylinders  or  spheres  at  steady  state  - constant  surface 
temperature  conditions]  can  be  transposed  logarithmically  to 


TC 


E 


R In  fPAlqZE 

IaRT26 


(6) 


and  the  critical  temperature  T~  can  be  iteratively  calculated  if  the 
physical  (p,  a,  6)  and  chemical  (Q,  Z,  E,  \)  parameters  are  known. 

Rogers  has  neatly  demonstrated9  the  general  validity  of  Eq.  (6)  by 
calculating  Tc  from  differential  scanning  calorimetric  determinations 
of  chemical  parameters  and  comparing  these  results  with  an  independ- 
ently measured  Tc.  The  excellent  agreement  between  experimental  and 
calculated  Tc  values  is  reproduced  in  Table  II. 


12Shouman,  A.R.,  Donaldson,  A.B.  & Tsao,  H.Y.,  Exact  Solution  to  the 
One  Dimensional  Stationary  Energy  Equation  for  a Self-Heating  Slab, 
Combustion  and  Flame  23  (1974)  17-28. 

90p  ait. 
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TABLE  II.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  CALCULATED 
CRITICAL  TEMPERATURES  (Tm) 


Tm  (°C) 

Values  used 

a 

P 

Q 

Z 

E 

AxlO4 

Exp 

Calc 

(cm) 

(g/cc) 

(cal/g) 

(sec-1 ) 

(kcal/mol ) 

(cal/cm°Csec) 

HMX 

253-255 

253 

0.033 

1 .81 

500 

5*1019 

52.7 

7.0 

RDX 

215-217 

217 

0.035 

1.72 

500 

2.02*1018 

47.1 

2.5 

TNT 

287-289 

291 

0.038 

1 .57 

300 

2.51 xlO11 

34.4 

5.0 

PETN 

200-203 

196 

0.034 

1.74 

300 

6. 3*1 019 

47.0 

6.0 

TATB 

331-332 

334 

0.033 

1.84 

600 

3.18xl019 

59.9 

10.0 

DATB 

320-323 

323 

0.035 

1 .74 

300 

1 .17x1 015 

46.3 

6.0 

BTF 

248-251 

275 

0.033 

1 .81 

600 

4.11 xlO12 

37.2 

5.0 

NQ 

200-204 

204 

0.039 

1 .63 

500 

2.84xl07 

20.9 

5.0 

PATO 

280-282 

288 

0.037 

1.70 

500 

1 .51xl010 

32.2 

3.0 

HNS 

320-321 

316 

0.037 

1 .65 

500 

1.53xl09 

30.3 

5.0 

The  important  point  to  be  learned  from  Roger's  work  is  that  good 
kinetic  constant  data  is  attainable  at  or  near  the  critical  temperature 
in  many  cases  and,  consequently,  Tc  values  can  be  calculated  with  some 
confidence.  But,  as  will  be  discussed  later  in  Section  IV,  unfortu- 
nately, this  is  not  necessarily  true  of  time-to-explosion  calculations. 

C.  Time-to-Explosion 

Again  referring  to  Eq.  (1),  in  those  special  cases  where 
there  is  no  transfer  of  heat  from  the  mass  (adiabatic  conditions), 
the  term  xv2T  approaches  zero  and  Eq.  (1)  reduces  to 


■ 13 

I 


NWSY  TR  76-2 


= QZ£-E/RT 


(7) 


Also,  again  making  the  same  simplifying  assumptions  used  for  deriving 
Eq.  (4),  namely  fy  « r|^|i  - and  e=  ^(T-To)  , Eq.  (7)  can  be 
expressed  as 


te 


cRTo„  E/RT 
QZE  * 


0e 


fc~Gde 


(8) 


Since  it  can  be  shown  that  e * 1 , the  time  to  ignition  or  explosion, 
te,  is  then  approximated  to  be 


te  - Si^E/KTo 
e QZE 


(9) 


where  T0  is  the  initial  temperature  of  the  adiabatical ly  controlled 
system. 

Eq.  (7)  is  often  conveniently  used  in  the  form  In  ^ = In  - E/RT 

which  is  analogous  to  the  empirical  equation  In  te  = A + B/T  that 
fits  adiabatic  temperature- time  profiles  for  many  explosives.  However, 
many  other  systems  do  not  fit  this  simple  linear  relationship. 
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Prior  to  the  work  of  Shouman,  et  al . ,12  Zinn  and  Rogers13  had  shown 
that  from  plots  of  the  nonlinear  relationship 


log10  WT  = function 


(10) 


times-to-explosion  values,  te,  could  be  read  from  the  curve  for 
specific  sample  surface  temperatures,  Ts,  where  E*  is  the  activation 
energy  of  the  rate  controlling  reaction  and  t = ~y~C.  Anderson14 
approximated  the  function  represented  by  Eq.  (10)  with  two  empirical 
straight  line  equations  limited  to  different  ranges  of  the  function, 
i.e. , 


log10  te/t  = 0.1400  - 0.0875  E* 


1_ 

Tc 


(ID 


for  te/r  in  the  range  0.17  to  1.2  and 


log10  te/t  = 1.134  - 0.1824  E*|^  - (12) 

for  te/x  in  the  range  0.001  to  0.17. 

Figure  4 reproduces  Anderson's  fit  of  his  empirical  equations  to  the 
Zinn/Rogers'  nonlinear  curve. 

120p  ait. 

13Zinn,  J.  & Rogers,  R.N.,  J.  Phys.  Chem.  66  (1962)  2646. 
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Table  III  is  constructed  from  data  presented  in  three  tables  by 
Sawyer  and  Williams5  who  calculated  critical  temperatures  and  times- 
to-explosion  for  the  ordnance  items  listed  using  Eq's.  (6),  (11)  and 
(12),  and  physical  and  kinetic  parameters  cited  by  Anderson14  and 
Joyner.15  Appropriate  6cr  values  are  used  for  those  charges  that 
were  not  simple,  solid  slabs,  cylinders  or  spheres.  Thus,  as  an  ex- 
ample, for  hollow  cylindrical  charges,  a corresponding  5cr  was  used 
in  the  calculations. 

The  5*10~9  probabilities  were  derived  statistically  by  Sawyer  and 
Williams  from  tables16  and  the  equation 


, _ 1 -SM 

A " SM(CV)  ' 


1 .84 


(13) 


where  Z = standard  normal  deviate 

CV  = coefficient  of  variation  in  thermal  response  data 
SM  * time-to-explosi on/exposure  time. 


Anderson's  equations  had  been  verified  experimentally  in  several 
cases  and,  on  the  face  of  it,  would  seem  to  imply  that  the  same 
kinetic  parameters  apply  in  Eq.  (6)  as  in  Eq's.  (11)  and  (12)  and 
that,  therefore,  one  might  expect  a similar  degree  of  confidence  in 
calculated  times-to-explosion  values  te,  as  in  those  for  critical 
temperature,  Tc.  That  this  is  often  far  from  the  case  comprises  the 
remainder  of  this  report. 

IV.  APPLICATIONS  AND  LIMITATIONS  OF  THERMOCHEMICAL  EXPLOSION  THEORY 
A.  Evaluation  of  ABL  Data 


Table  III  suffers  from  a serious  flaw  due  to  simplifying 
assumptions  in  its  construction.  The  originators  state  that 
"...quantitative  engineering  analysis  was  made  to  determine. . .exposure 
times/temperatures  corresponding  to  an  acceptable  risk  level."  This 
statement  does  not  heed  the  warnings  in  some  of  their  prime  data 
sources. 


5Op  oit. 

14Ariderson,  C.M.,  NWC  TP  5245,  Cookoff  Prediction  of  Explosives  by 
Thermal  Analysis  Techniques,  Jul  1972. 

15Jo.yner,  T.B.,  NWC  TP  4709,  Thermal  Decomposition  of  Explosives  (U), 
Parts  1,2,4&9,  dtd  Mar  1969,  Apr  1969,  Jun  1970  & Oct  1973  (resp) 
CONFIDENTIAL. 

16J.  Am.  Statistical  Assn.,  Jun  1970,  p 637. 
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TABLE  III.  ABL  CALCULATED  TIMES-TO-EXPLOSION  FOR  ORDNANCE  ITEMS  AT  VARIOUS  TEMPERATURES 


Critical 

Contents 

Times-to-explosiona  (hrs) 

Item 

temp 

(°F) 

(molten) 

21 2°F 

228°F 

241  °F 

250 

°F 

259 

°F 

A 

B 

A 

B 

A 

B 

A 

B 

A 

B 

Mines 
Mk  25 

221 

HBX-1 /Asphalt 

>1 26 . 0b 

>81 .0 

114.0 

74.0 

74.6 

49.0 

55.8 

37.0 

42.0 

28.0 

Mk  25 

180 

TNT/Asphalt 

127.0 

82.0 

84.8 

55.0 

62.6 

41 .0 

50.7 

33.0 

41 .3 

27.0 

Mk  39 

221 

HBX-1 /Asphalt 

>126.0 

>81 .0 

114.0 

73.0 

74.6 

48.0 

55.8 

36.0 

42.0 

27.0 

Mk  50 

239 

HBX-3/Asphal t 

>29.0 

>18.0 

>29.0 

>18.0 

>29.0 

>18.0 

23.5 

15.0 

17.9 

11.5 

Torp  WH 

Mk  17 

221 

HBX-1 /Asphalt 

>111 .0 

>72.0 

99.8 

65.0 

65.6 

43.0 

49.0 

31 .0 

36.9 

24.0 

Mk  37 

244 

HBX-3/Asphalt 

>19.2 

>12.5 

>19.2 

>12.5 

>19.2 

>12.5 

18.6 

12.0 

14.0 

9.0 

Dep  Chg 

Mk  8 

205 

TNT/ Asphalt 

157.0 

100.0 

104.0 

67.0 

77.0 

50.0 

62.4 

41.0 

50.8 

33.0 

Bombs 
Mk  82 

239 

H-6/Asphalt 

>26.3 

>17.0 

>26.3 

>17.0 

>26.3 

>17.0 

21.3 

10.5 

16.1 

5.7 

Mk  82 

203 

TNT/ Asphalt 

27.1 

17.5 

18.1 

11.8 

13.3 

8.7 

10.8 

7.0 

8.8 

5.7 

A = 0.5  probability  at  50%  confidence. 

B = 5x10” 9 probability  at  95%  confidence. 

aTime  zero  is  when  item  is  first  subjected  to  the  given  temperature. 

^Anderson's  technique  of  calculating  times-to-explosion  requires  the  use  of  certain  parameters 
which  are  extracted  from  experimental  data.  For  several  items,  times-to-explosion  could  not 
be  determined  at  certain  temperatures  because  required  parameters  were  outside  the  given  range. 
For  such  cases,  times-to-explosion  were  calculated  using  end-of-the-range  parameters.  The 
symbol  > appears  before  each  number  so  calculated,  indicating  that  the  time-to-explosion  is 
actually  greater  for  the  given  temperature  than  indicated  by  that  number  given. 
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The  questionable  assumptions  were: 

• That  kinetic  parameters  derived  from  "chemically  pure" 
systems  would  be  the  same  for  the  aged  bomb  compositions. 

• That  the  kinetic  parameters  derived  for  the  critical 
temperature  calculations  are  the  same  rate  controlling 
factors  as  those  for  times-to-explosion. 

• That  no  side,  autocatalytic  or  other  competing  reactions 
are  occurring  in  any  of  the  systems  treated. 

In  short,  our  objective  here  is  to  reemphasize  that  a better  under- 
standing of  the  chemical  ambiguities  involved  in  explosive  systems 
is  essential  to  meaningful  data  reduction  by  mathematicians,  stat- 
isticians and  heat  transfer  specialists,  as  well  as  by  explosive 
design  and  process  engineers.  Section  IV. B summarizes  relevant 
caveats  that  apparently  were  not  taken  into  account. 

B . Evaluation  of  the  Thermal  Explosion  Literature 

1 . Autocatalysis 

Longwell6  warns  of  the  autocatalytic  decomposition  of 
nitrocellulose  (used  in  D-2  waxes  and  all  HBX-type  explosive  compo- 
sitions) and  points  out  that  calculated  values  of  Tg  and  te  "...should 
be  used  with  margins  of  safety  dictated  by  engineering  judgments  that 
are  strongly  influenced  by  the  magnitude  of  the  cost  of  being  wrong." 

2.  Solid-state  decomposition 


Rogers* 9 *  emphasizes  the  key  points:  1)  that  critical 

temperatures  are  determined  by  the  maximum  rate  attainable  in  the 
condensed  phase  reaction  for  the  explosive  in  question,  and  2)  that 
kinetic  constants  from  these  rapid  reactions  can  be  measured  "quite 
accurately"  but  that  times-to-explosion  may  be  largely  determined  by 
difficult-to-measure,  premolten,  solid-state  reactions  which  "...can 
vary  tremendously  with  changes  in  purity  and  crystal  perfection." 

3.  Catalysis 

Joyner15  speaks  of  the  catalytic  role  of  iron  oxides 
(rust)  causing  as  much  as  tenfold  acceleration  rates  in  TNT-based 
systems  and  extrapolation  of  his  TNT/asphalt  (clean)  data  shows  that 
at  120°C  there  is  already  a 1 percent  decomposition  in  only  3 hours. 


6<9p  cit. 

9 Op  cit. 

15 Op  cit. 
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He  also  warns  of  the  many  potential  side  reactions  of  past  and 
current  weapons  systems  containing  some  or  all  of  the  following 
ingredients  in  various  states  of  purity:  TNT,  RDX,  hot  melt,  wax, 

CaCl2»  H2O,  A1 , rust,  NH4NO3,  KNO3,  paint,  etc. 

Joyner's17  basic  warning  is  that  in  order  to  substantiate  theory, 
most  experimental  data  is  generated  from  highly  purified  materials 
where  reaction  products  are  usually  free  to  escape.  Until  fairly 
recently,  no  serious  efforts  have  been  made  to  simulate  "real  world" 
reaction  conditions.  The  kinetic  constants  generated  by  Joyner,  and 
used  by  Anderson  and  ABL,  were  those  derived  from  highly  purified 
components  with  escaping  reaction  products. 

Urbanski18  also  describes  accelerated  reaction  rates  for  TNT-metal 
systems. 


4.  Temperature  control 

The  AMCP  handbook7  warns  of  the  importance  of  both 
accurate  temperatures  and  close  temperature  controls  for  any  heating 
operations  such  as  melting,  drying,  etc.  Even  assuming  that  none  of 
the  problems  cited  in  B.l,  2 and  3 exist  in  aged  ordnance,  the  sta- 
tistical calculations  for  5xl0~9  probabilities  of  explosion  made  by 
ABL  give  the  explosive  process  engineer  a false  sense  of  security. 
Figure  5,  generated  from  their  calculated  values  in  Table  III  for  the 
TNT-asphalt  Mk  25  Mine,  illustrates  this  point.  It  can  be  seen  that 
relatively  small  temperature  changes,  i.e.,  only  about  10°C,  can 
lead  to  a 50%  probability  of  explosion  rather  than  a 5x1 Q~9  probability 
for  the  same  systems  heated  for  the  same  time  period. 

That  such  temperature  fluctuations  are  far  from  an  unlikely  occurrence 
in  production  operations  is  discussed  next. 

5.  Explosion  in  the  TNT  melt  shop  at  Muiden,  Holland 

Groothuizen19  and  his  coworkers  have  reported  on  their 
investigation  of  the  KNSF  melt  shop  explosion  in  June  1966.  While 
the  TNT  reclamation  process  going  on  prior  to  the  explosion  does  not 
completely  parallel  the  autoclave  melt-out  process,  there  are,  none- 
theless, uncomfortable  similarities. 


17Joyner,  T.B.,  private  communication,  Dec  1975. 

18Urbanski , T.,  Chemistry  & Technology  of  Explosives,  Vol  I,  Pergamon 
Press  Ltd. , 1967. 

7Op  cit. 

19Groothuizen , T.M.,  Lindeijer,  E.W.  & Pasman,  H.J.,  Explosivestoffe 
5 (1970)  97. 
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In  the  KNSF  melt  shop,  the  melt  tank  was  an  agitated,  steam-jacketed 
vessel  with  bottom  and  side  drain  cocks  and  an  inside  diameter  of 
3 feet.  The  melt  tank  was  used  to  reclaim  TNT  from  various  sources 
for  reloading  into  other  weapons  and  contained  about  1000  pounds  of 
TNT  at  the  time  the  explosion  occurred. 

Unfortunately,  the  operating  conditions  prior  to  the  explosion  were 
not  completely  known.  For  examples:  the  explosion  occurred  at  4:00 

p.m.  on  2 June  with  the  melt  tank  being  put  in  operation  on  31  May 
but  it  is  not  clear  whether  it  was  continuously  being  heated  from 
31  May  to  2 June  or  only  on  one  shift  per  day,  temperatures  may  have 
gone  as  high  as  150°C,  or  even  higher,  for  reasons  not  clearly  listed; 
various  contaminants  such  as  cardboard,  KC1 , NH4NO3,  KNO3,  etc.,  may 
or  may  not  have  been  present  in  the  scrap  being  melted;  the  agitator 
was  not  in  operation  at  explosion  time. 

Nonetheless,  what  is  important  were  the  results  of  the  post-explosion 
investigation.  Four  samples  were  formulated  from  TNT  and  supposed 
contaminants  and  held  adiabatically  while  the  heat  evolution  of 
reaction  was  measured  against  time.  Samples  were: 

(1)  Pure  TNT 

(2)  TNT  + 2%  NH4NO3  + 1%  cardboard  + trace  KC1 

(3)  TNT  + 0.2%  KNO3 

(4)  TNT  + 0.5%  propellant  + 0.5%  soap  powder  + 0.5%  carboard 
+ 0.5%  Fe  + 2%  NH4NO3  + trace  KC1 

For  the  first  three  cases,  "...heat  generated  increased  significantly 

with  temperature  (2.3  to  3 times  per  10°C) At  the  same  temperature 

(180°C)  the  respective  heat  generations  for  samples  (2)  and  (3)  were 
20  and  35  times  higher  than  the  pure  TNT  (1)."  In  sample  (4),  which 
in  many  ways  most  resembles  typical  aged  Navy  ordnance  that  will  be 
demilitarized  in  the  future,  only  one  measuring  point  was  obtained. 
"The  reason  was  that  5 hours  after  the  test  started,  the  mixture 
suddenly  ignited.  Prior  to  that,  the  temperature  had  risen  only  a 
few  degrees.  The  temperature  just  before  ignition  was  134°  to  135°C." 

Other  findings  by  Groothuizen,  ct  at.,  are  instructive.  In  their 
investigations  of  cardboard  contaminants,  they  concluded  that  "...it 
is  indeed  possible  that  by  absorption  of  TNT  on  a solid  having  a large 
specific  surface,  the  rate  of  decomposition  is  greatly  increased. 

This  appeared  from  tests  with  active  carbon  as  a contaminant  in  which 
the  TNT  ignited  15  minutes  after  the  mixture  had  reached  the  oven 
temperature  of  135°C." 
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They  hypothesized  that  the  abrupt  ignition  of  TNT  when  cardboard  is 
present  is  due  to  the  formation  of  a solid  mass  of  TNT  together  with 
an  insoluble  decomposition  product.  In  this  mass,  convection  is 
impossible  and  the  slow  dissipation  of  heat  by  conduction  only  can 
lead  to  local  self-heating  and  hot  spots.  This  explanation  is  sup- 
ported by  the  fact  that  abrupt  ignition  also  took  place  when  glass 
wool  was  substituted  for  cardboard. 

Also,  in  an  earlier  study,  this  same  laboratory20  had  made  another 
highly  significant  finding  - the  effect  of  ageing. 

In  brief,  measurements  were  made  of  the  adiabatic  self-heating  of 
TNT  versus  time  for  both  pure  TNT  and  for  TNT  that  had  been  aged  by 
successive  melting,  cooling  and  remelting,  etc.  Data  obtained 
fitted  the  familiar  Arrhenius  type  equation: 


q = Cat“Ea/pT 


(14) 


where  q = heat  generated  watts/kg 
R = gas  constant  cal/mol°K 
Ea  = apparent  energy  of  activation  cal/mol 
Ca  = pre-exponential  factor  watts/kg 
T = temperature  °K 


The  data  obtained  is  highly  instructive  and  is  reproduced  in  Table  IV. 


TABLE  IV.  ADIABATIC  SELF-HEATING  OF  TNT 


Sample 

Temp  range 
(°C) 

Depth  of  con- 
version as  °C 
self-heating 

LogioCa 

(w/kg) 

Ea 

(Kcal/mol ) 

Pure  TNT 

145-196 

0-51 

20.60 

43.0 

Remelted  TNT 

150-160 

51-61 

15.46 

31.3 

Re-remelted  TNT 

126-144 

61-79 

13.51 

27.2 

20Pasman,  H.J.,  Groothuizen,  T.M.  & Vermeulen,  C.M.,  Explosivestoffe 
7 (1969). 
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The  striking  increase  in  depth  of  conversion,  61-79  at  126°  to 
144°C  for  re-remelted  TNT  versus  0-51  at  145°  to  195°C  for  pure  TNT, 
coupled  with  the  marked  decrease  in  the  exponential  apparent  energy 
of  activation,  Ea,  raises  serious  questions  about  the  advisability 
of  melt-recast  processes  without  some  intermediate  purification  of 
the  explosive  system.  This  is  in  keeping  with  the  tenfold  decrease 
in  vacuum  stability  and  the  methyl  violet  reaction  previously  de- 
scribed in  the  NAVWPNSTA  Yorktown  work  in  Section  II. 

In  addition,  Groothuizen,  et  at.,  observed  gassing  prior  to  ignition 
or  explosion  in  their  controlled  experiments.  This  factor  also  par- 
allels the  NAVWPNSTA  Yorktown  experiment.  When  the  original  auto- 
clave experiment  was  run  at  NAD  Hawthorne,  gassing  was  also  noted 
but  was  attributed  to  trapped  air. 

As  a final  note  of  warning,  it  is  a current  Army  practice  to  subject 
major  caliber  projectiles  to  long  (several  days)  heat  soaking  at 
elevated  temperatures  to  heal  cracks  and  case  separations  of  their 
explosive  charges. 

Potential  decomposition  problems  are  obviously  present  and  it  is 
recommended  that  further  studies  along  the  lines  of  Groothuizen  and 
his  co-workers  be  carried  out  to  assure  that  acceptable  safety  margins 
exist. 

V.  SUMMARY 

• During  the  autoclave  melt-out  of  aged  Mk  25  Mines,  decomposi- 
tion reactions  began  to  take  place  at  times  well  below  those 
predicted  from  thermochemical  theory  for  pure  components  in 
the  system. 

• Ample  warnings  exist  in  the  thermochemistry  of  explosion 
literature  testifying  to  the  fact  that  many  side  and  solid- 
state  reaction  kinetics  may  be  the  rate  determining  steps 
for  times-to-explosion  in  less  than  pure  systems  of  the 
above  types. 

• A purely  statistical  estimate  of  probabilities  of  explosion 
even  for  pure  systems  can  be  beguiling  without  mention  of 
the  tremendous  impact  of  small  changes  in  temperature  on 
overall  explosion  times. 

• The  differences  between  theoretically  calculated  and  experi- 
mentally determined  values  for  critical  temperatures  and 
times-to-explosion  have  been  steadily  narrowed  for  pure 
systems  by  a succession  of  outstanding  investigations. 

However,  a better  understanding  and  more  kinetic  data  are 
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required  of  impure  reaction  systems  before  any  quantitative 
thermal  explosion  predictions  can  be  made  about  them.  This 
is  particularly  true  when  prolonged  and/or  repeated  heating 
of  any  kind  is  involved. 
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